Unified picture of superfluidity: Prom Bogoliubov's approximation to Popov's 

hydrodynamic theory 
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Using a non-perturbative renormalization-group technique, we compute the momentum and fre- 
quency dependence of the anomalous self-energy and the one-particle spectral function of two- 
dimensional interacting bosons at zero temperature. Below a characteristic momentum scale kc, 
where the Bogoliubov approximation breaks down, the anomalous self-energy develops a square root 
singularity and the Goldstone mode of the superfluid phase (Bogoliubov sound mode) coexists with 
a continuum of excitations, in agreement with the predictions of Popov's hydrodynamic theory. 
Thus our results provide a unified picture of superfluidity in interacting boson systems and connect 
Bogoliubov's theory (valid for momenta larger than kc) to Popov's hydrodynamic approach. 

PACS numbers: 05.30.Jp,03.75.Kk,05.10.Cc 



Following the success of the Bogoliubov theory [ij 
in providing a microscopic explanation of superfluidity, 
much theoretical work has been devoted to the calcula- 
tion of the one-particle Green function of Bose super- 
fluids 0- Early attempts to improve the Bogoliubov 
approximation however encountered difficulties due to 
a singular perturbation theory plagued by infrared di- 
vergences [3[. In the 70s, Nepomnyashchii and Nepom- 
nyashchii proved that the anomalous self-energy (the 
main quantity determining the one-particle Green func- 
tion) vanishes at zero frequency and momentum in di- 
mension d < 3 [1,0]. This exact result shows that the Bo- 
goliubov approximation, where the linear spectrum and 
the superfluidity rely on a flnite value of the anomalous 
self-energy, breaks down at low energy. The vanishing 
of the anomalous self-energy in the infrared limit has a 
definite physical origin; it is due to the coupling between 
transverse (phase) and longitudinal fluctuations and the 
resulting divergence of the longitudinal susceptibility - a 
general phenomenon in systems with a continuous broken 
symmetry Q. 

An alternative approach to superfluidity, based on a 
phase-amplitude representation of the boson fleld, has 
been proposed by Popov 01 ■ This approach is free of in- 
frared singularity, but restricted to the (low-momentum) 
hydrodynamic regime and therefore does not allow to 
study the higher-momentum regime where the Bogoli- 
ubov approximation is valid. Nevertheless, Popov's the- 
ory agrees with the exact result and the asymptotic 
low-energy behavior obtained by Nepomnyashchii and 
co-workers [1, [^. Furthermore, the expression of the 
anomalous self-energy obtained by Nepomnyashchii et al. 
and Popov in the low-energy limit yields a continuum of 
(one-particle) excitations coexisting with the Bogoliubov 
sound mode [lO| , in marked contrast with the Bogoliubov 
theory where the sound mode is the sole excitation at low 
energies. 



Although the breakdown of the Bogoliubov approxi- 
mation in d < 3 is now well understood within the renor- 
malization group approach 11, l^, 13, 14 1, no theoretical 
framework has given a unifled description, from high to 
low energies, of superfluidity in interacting boson sys- 
tems. Taking advantage of recent progress in the non- 
perturbative renormalization group (NPRG), we have 
calculated the momentum and frequency dependence of 
the anomalous self-energy of interacting bosons. Our re- 
sults provide a unifled description of superfluidity that 
encompasses both Bogoliubov's theory and Popov's hy- 
drodynamic approach. 

In this Letter, wc focus on a two-dimensional sys- 
tem [l^. We show that the Bogoliubov approxima- 
tion breaks down at a characteristic momentum scale 
kc which, for weak boson-boson interactions, is much 
smaller than the inverse healing length kh [iBl- Moreover, 
the anomalous self-energy Ean(p, 'j-') becomes singular be- 
low kc and induces a continuum of excitations which co- 
exists with the Bogoliubov sound mode, in agreement 
with the predictions of Popov's hydrodynamic approach. 
We compute the longitudinal (one-particle) spectral func- 
tion A\i{p,Lu) and discuss its dependence on |p|/fcG [l7| . 

We consider two-dimensional interacting bosons at 
zero temperature, with the action 



S = dx 



{h = kB = 1 throughout 
is a bosonic (complex) fleld. 



the Letter), where ip{x) 
X = (r, t), and J dx = 
JI^ dr J d^r. r S [0, /?] is an imaginary time, /3 — s- oo 
the inverse temperature, and denotes the chemical po- 
tential. The interaction is assumed to be local in space 
and the model is regularized by a momentum cutoff A. It 
is convenient to write the boson fleld as i/; = ■^{i'l+i^^) 
with and ■02 real. We assume the dimensionless cou- 



2 



pling parameter 2gm to be much smaller than unity 
(weak coupling limit). 

The excitation spectrum can be obtained from the one- 
particle Green function 



'-G,,{p-n)+\5, 



eijG\t{p; n) 



Gtt(p; n) 



(2) 



or its inverse, the two-point vertex 

rif (p; 0) = k,j^A{p] n) + (j),<j)jTB{p] n) + eyTc(p; n), 

(3) 

where is the antisymmetric tensor and p — (p, iuo) 
(with uj a Matsubara frequency), (pi = (ipi^x)) (i = 
1,2) is assumed to be space and time independent, and 
n = (0f -I- 02)/2 denotes the condensate density. Tak- 
ing advantage of the (global) gauge invariance of the ac- 
tion (P), we have introduced different Green functions 
for transverse (phase) and longitudinal fluctuations wrt 
the constant field (01 ,4>2)- Formally, the two-point ver- 

(2) 

tex r^^- ' {p; (j)) can be defined as the second-order func- 
tional derivative of the effective action r[c/)] (the generat- 
ing functional of one-particle irreducible vertices). The 
NPRG procedure is set up by adding to the action ([T|) 
an infrared regulator AS"/,. = X]p V'*(p)-Rfc(p)^(p) which 
suppresses fluctuations with momenta below a charac- 
teristic scale k but leaves the high-momenta modes unaf- 
fected (for a review, see Ref. [ill)- The dependence of the 
effective action on k is given by Wetterich's equation [l^ 

dkVj,[4>] = ^Tr {5fci?fc(rf [0] + RkY'] . (4) 

In Fourier space, the trace involves a sum over frequencies 
and momenta as well as a trace over the two components 
of the field cj) = {(j)i,(j)2)- For ~ A, the regulator ASk[(l)] 
suppresses fluctuations and the mean-fleld theory, where 
the effective action TA[(j)] reduces to the microscopic ac- 
tion S[(p], becomes exact thus reproducing the results of 
the Bogoliubov approximation. On the other hand for 
/c = 0, provided that Rk{p) vanishes, Tk[(j)] gives the ef- 
fective action of the original model H]) from which we 
can deduce F^^^ and the one-particle Green function. 

We make the following two simplifications: i) we ne- 
glect the field dependence of Ta^k{p',n) {a = A,B,C) 
which we approximate by its value at the actual (fc- 
dependent) condensate density no,fc- '^o,fc is obtained 
from the minimum of the effective potential Uk{n) = 
{f3V)-^rk[^] with (/) = (72^,0) {V is the volume of the 
system), ii) We approximate the effective potential by 



Ukin) = C/fc(rio,fc) + ^{n ~ no,k) 



(5) 



where Aj^^A = g. Our approach follows the NPRG 
scheme recently proposed by Blaizot, Mcndcz-Galain and 
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FIG. 1: (Color online) \/g, Zc, Va/VX and c/c|fe^o vs. 
\n{kG/k), where kc = {gmy^n / An , for n = 0.01 and 
2mg = 0.1 {ln{kG /kh) — —5.87). The inset shows kc vs. 2mg 
obtained from the criterion VaUg ~ ^1/2 (the green solid line 
IS a fittoifcG tx {2mgf''^). AUfi sures are obtained with A = 1, 
2m = 1 and the regulator i?fc(p) = ZA,kep/{e^ - 1). 



Wschebor and others [20|, |2l|, |22| with a truncation in 
fields to lowest non-trivial order [23| . 

Previous NPRG studies of interacting bosons assumed 
a simple form of the effective action F;; [0] with local and 



0{V\dr,dl) terms [12|, llj, Uj . In our formalism, this 
amounts to expanding the vertices V a , F_b , Fc in powers 
of p and 

Fa(p) = VAW^ + Z^cp, Fb(p) = A, Vc{p) = Zc^ (6) 

(we drop the k index to alleviate the notation), where 
Ep = 12m, is the dispersion of the free bosons. The 
initial conditions (fc = A), Za = Zc ^ I, Va = 0, ^ ~ g 
and no = reproduce the Green function Gij{p) of 
the Bogoliubov approximation. Equations (|2l6p yield a 
low-energy mode uj ~ c|p|, with (fc-dependent) velocity 



ZA/{2m) 
VA + Z^/{2Xno)J ' 



(7) 



and a superfluid density n,, = Z^ng = n where n is 
the mean boson density [ij, In the weak-coupling 
limit, uq and Za are weakly renormalized and Uq ~ ^J'/g, 
Z\ ~ 1 (with TT-Q = no|/£=Oj etc.). On the other hand 
Zc ~ fc and A ~ fc vanish when fc ^ while Va takes 
a finite value VI (Fig. [1]) [H, [i3| . The anomalous self- 
energy San(0,0) = noFs(0, 0) = noA therefore vanishes 
for fc ^ in agreement with the exact result l^. The 
existence of a linear spectrum is then due to the rela- 
tivistic form of the action which emerges at low energy 



and ^ > 0) [la, lij. The character- 



istic momentum scale kc ("Ginzburg" scale) at which 
the Bogoliubov approximation breaks down can be de- 
fined by the criterion VaUg = ^a/^- In the weak cou- 
pling limit 2gm ^ 1, it is found to be proportional to 
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FIG. 2: (Color online) Real and imaginary parts of the re- 
tarded vertex r§(p,aj) = T,^^{p,ijj)/no for various values 
of IpI ranging from O.Skc up to ITOkc ~ kh/2 (n — 0.01, 
2mg — 0.1 and fc = 0). The Bogoliubov approximation corre- 
sponds to r^{p,iu) = 5 = 0.1. (3?[ri] > and 9[r||] < 0.) 
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FIG. 3: (Color online) Same as figure [2] but for |p|/fcG — 
0.036. The dotted lines show the analytical result (|9} obtained 
from the approximation ([8]). The latter is shown by the blue 
dashed-dotted fine while the blue crosses give the numerical 
results for rs(p,i(j;) vs. the rescaled Matsubara frequency 
[ij/(c*|p|). The inset shows r§(p,a;) on a much larger energy 
scale {uj < 10c* fee ~ 280c* |p|; the square root singularity at 
UJ <ti c*kG is not visible on this plot). 



-\/ {gm)^n ^ gmkh (fc/i = \J2mgn is the inverse heal- 
ing length below which the spectrum becomes linear) 
in agreement with a simple estimate based on the (one- 
loop) perturbative correction to the Bogoliubov approx- 
imation [ll[ (sec the inset in Fig.[T]). In practice, we use 
the definition kQ = ^ {gm)-^fi/ (An). 

Now we discuss the momentum and frequency depen- 
dence of the vertices beyond their derivative expansion 
(|6]). To simplify the numerical evaluation of dk^a,k{p) 
(a = A,B,C), we approximate the propagators enter- 
ing the flow equations using This type of approx- 
imation has been shown to be very accurate for classi- 
cal systems 2J]. It is sufficient to obtain the anoma- 
lous self-energy San(p) = noTBip) but is less rehable 



for the calculation of the (small) damping terms aris- 
ing from Ta and Tc [l^. In practice, wc compute the 
vertices Ta{p,iuj) for typically 100 frequency points and 
then use a Pade approximant to obtain the retarded 
part r^(p,w) = Taip,uj + iO+) While the Bogoli- 

ubov result r^(p,tt') = g is a good approximation when 
IpI » kc, r|^(p,Lj) develops a strong frequency depen- 
dence for IpI < kc (Fig. [2|). In the Hmit |p| <^ kc and 
|a;| <C c*kG, the vertices Ta and Tc are very well ap- 
proximated by their low-energy limit ([6|), while Tb{p) is 
well fitted by a square-root singularity. 



TB{p,iuj)o^C^yu;^ + {c*p) 



(8) 



with C a p-independent constant (the Bogoliubov result 
Tg{p,uj) ~ 5 is nevertheless reproduced at high energies 
\uj\ > c*kc) (Fig. [3]). For the retarded part, equation ([8]) 
implies 



ri(p,^) 



Ce{c*\p\-to)^ic*py-u^ 
-iC6i(cj-c*|p|)Vw2- (c*p)2 (9) 



{9{x) is the step function) when |p|,w/c* <C kc (we dis- 
cuss only oj > 0). The square root singularity ^ of 
the anomalous self-energy noTg{p,uj) agrees with the 
result obtained from diagrammatic resummations Q or 
the predictions of the hydrodynamic approach Q in the 
limit {p,oj) 0. As shown in Fig. [3l this singularity 
is very well reproduced by the result deduced from the 
Fade approximant. Thus our results interpolate between 
Bogohubov's approximation (|p| kc) and Popov's hy- 
drodynamic theory (|p| ^ kc). 

In the low-energy limit p,cj — + 0, Tg{p,uj) is of order 
IpI, |a;|, while r^(p, uj) and r^(p, uj)'^ are of order p^, uj'^. 
The onc-particle Green function then becomes 



G^^(p,w) 
Gi{p,Lu) 

G5(p,c.) 



T'^ip,u;) + 2n*Tlip,LU 



D^{p,u;) 



riiP,u;y 



D^(p,tj) 

R( 



2nSr«(p,u;)' 



(10) 



r^(p,c.) 



D'^{p,uj)' 



where D^ip,u;) = r2(p, c^)[r«(p, c^) + 2nSrf (p,u;)] + 
rg(p,w)2. From jg and one concludes that the 

transverse spectral function. 



^tt(p,w) 



-3[Gt1(p,c.) 



S{cO'C*\p\ 
2FJ:c*|p| 



(11) 



(for UJ > 0), exhibits a Dirac-likc peak at the Bogoli- 
ubov mode frequency to = c*|p|, in very good agreement 
with the result obtained from the Pade approximant [25| . 
From dH]) and pO|) , we deduce that the longitudinal spec- 
tral function 



^ii(p,^) 



1 



e{uj^c*\p\) 



2TTn*C ^Cj2 - (c*p)2 



(12) 
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FIG. 4: (Color online) Longitudinal spectral function 
Aii{p,Lj) for |p|/A;g — 0.036 and the same parameters as in 
figure [3l The red solid line is the result obtained from the 
Fade approximant, while the green dashed line is obtained 
from the analytic expression (|12|l . The inset shows the ratio 
between Aii(p, uj) and the approximate result (|12|l on a larger 
energy scale. 



exhibits a continuum of excitations with a singularity 
at the Bogoliubov mode frequency uj = c*|p|, again in 
agreement with the prediction of the hydrodynamic ap- 
proach [l^. The analytic expression (|12p gives a good 
approximation to the result obtained from the Pade ap- 
proximant when IpI ^ kc and \uj\ < c*kc (Fig.H]). This 
defines the domain of vahdity of the Popov hydrody- 
namic approach. For |p| ^ kc, the continuum of excita- 
tions is suppressed and the longitudinal spectral function 
Aii{p,uj) reduces to a Dirac-like peak as in the Bogoli- 
ubov theory. 

It should be noticed that we have retained only the 
leading term Tb{p ^ 0) = X in the derivative expansion 
of Tsip) [Eq. This is made possible by the fact that 
r b{p = 0) = X = 0{k) is a. very large energy scale wrt 
Ta{p),^c{p) ^ foi' typical momentum and frequency 
|p|, a;/c = 0{k ). The (p,a;) dependence of Tb{p) A + 
C\JuP- + (cp)^ does not change the leading behavior of 
TbCp) = C'(fc) which acts as a large mass term in the 
vertices. It should also be noticed that the singularity of 
r^^) yields a similar singularity in higher-order vertices 0, 
i]: r(3) ~ .J^B and F^^) ^Vb- These singular terms 
are neglected in our approach and we have approximated 
F(3) ~ and F^'') A. Again this is justified by the 

fact that A = 0(fc). 

To conclude, we have obtained a unified description 
of superfluidity which is valid at all energy scales and 
connect Bogoliubov's theory to Popov's hydrodynamic 
approach. Our results reveal the fundamental role of 
the Ginzburg momentum scale kQ in interacting boson 
systems, kc sets the scale at which the Bogoliubov 
approximation breaks down and determines the region 
|p|, Iti'l/c* ^ kQ where the longitudinal spectral func- 
tion An(p,w) _ (c*p)2 [Eq. ([ni)], i.e. the do- 



main of validity of the hydrodynamic approach. From a 
more general perspective, our results also show that the 
NPRG is a very efficient tool to study strongly correlated 
quantum systems and in particular to compute spectral 
functions. 

The author would like to thank B. Delamotte and N. 
Wschebor for useful discussions. 
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